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I Abstract. We construct an explicit solution of the Cauchy initial value problem for the time- 

• dependent Schrodinger equation for a charged particle with a spin moving in a uniform magnetic field 

^ I and a perpendicular electric field varying with time. The corresponding Green function (propagator) 

is given in terms of elementary functions and certain integrals of the fields with a characteristic 
•/^ I function, which should be found as an analytic or numerical solution of the equation of motion for 

the classical oscillator with a time-dependent frequency. We discuss a particular solution of a related 
(~| ' nonlinear Schrodinger equation and some special and limiting cases are outlined. 

1. Introduction 
\0 • The time-dependent Schrodinger equation 

I: i^^-nm, (1.1) 

' can be solved using the time evolution operator given formally by 
O : U{t, to) = T (^exp j'^ H it') dt'^ ^ , (1.2) 



- '—I 

X 



where T is the time ordering operator which orders operators with larger times to the left [3] , [12] 
This unitary operator takes a state at time to to a state at time t, so that 



^P{x,t) = Uit,to)^ix,to). (1.3) 

The simplicity of these formulas is deceptive, since the time evolution operator can be found explic- 
itly as certain integral operator only in a few special cases. An important example of this source 
is the forced harmonic oscillator originally considered by Richard Feynman in his path integrals 
approach to the nonrelativistic quantum mechanics [7], [8], [9], [10], and [11]; see also [21j. Since 
then this problem and its special and limiting cases were discussed by many authors; see Ref. [1], 
[13], [16], [22], [21], [39] for the simple harmonic oscillator and Ref. [2], [6], [15], [27], [33] for the 
particle in a constant external field and references therein. It is worth noting that an exact solution 
of the n-dimensional time-dependent Schrodinger equation for certain modified oscillator is found 
in 
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In this Letter we construct the time evolution operator exphcitly in a general case of the one- 
dimensional Schrodinger equation when the Hamiltonian is an arbitrary quadratic form of the 
operator of coordinate and the operator of linear momentum; see equation (12. ip below. In this 
approach, all exactly solvable models, that have been cited above, are classified in terms of ele- 
mentary solutions of a characterization equation given by fl2.13p below. A particular solution of 
the corresponding nonlinear Schrodinger equation is obtained in a similar fashion. By separation 
of variables, we apply this method to another classical problem — the motion of a charged parti- 
cle with a spin in uniform magnetic and perpendicular electric fields that are changing with time. 
The corresponding Green function (or Feynman's propagator) is found as an elementary function of 
certain integrals of our characteristic function and the electromagnetic fields; see equation 06.141) be- 
low. Special cases of constant and hnear magnetic fields are discussed as examples. Moreover, these 
explicit solutions can also be useful when testing numerical methods of solving the time-dependent 
Schrodinger equation. 



2. Solution of a Cauchy Initial Value Problem 



The fundamental solution of the time-dependent Schrodinger equation with the quadratic Hamil- 
tonian of the form 

= -« (^) 1^ + ^ W ^'V^ - ' W ^^ + ^ W ^) - / (^) + '9 (t) ^, (2.1) 

where a{t) , b {t) , c{t) , d (t) , f {t) , and g {t) are given real- valued functions of time t only, can be 
found with the help of a familiar substitution 

^ = Ae'^ = A (t) e^^^^'J''*) (2.2) 

with 

A = A{t)= , ^. (2.3) 
A/27riyU (t) 

and 

S = S {x,y,t) = a (t) + (3 (t) + 7 (t) + 6 (t) x + e (t) y + k (t) , (2.4) 

where a{t) , (3 (t) , 7 (t) , 5 (t) , e (t) , and k, (t) are differentiable real-valued functions of time t only. 
Indeed, 

dS fdSV , . , .dS 



- = -a^-J -bx^ + fx + ig-cx)- (2.5) 
by choosing 

K. = a^ + d = 2a{t)a{t) + d{t). (2.6) 

Equating the coefficients of all admissible powers of x™?/" with < m + n < 2, gives the following 
system of ordinary differential equations 

dcy. 

— + b{t) + 2c{t)a + Aa{t)a'^ = 0, (2.7) 
d3 

-£ + {c{t)+4a{t)a{t))P = 0, (2.8) 
^ + a (t) (t) = 0, (2.9) 
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dS 

- + (c it) + 4a (t) a (t)) 6 = f{t) + 2a (t) g (t) , (2.10) 
at 

de 

- = {g{t)-2a{t)6{t))f3{t), (2.11) 

- = g{t)6{t)-a{t)6'{t), (2.12) 



where the first equation is the famihar Riccati nonhnear differential equation; see, for example, 
[32] . [lO] and references therein. Substitution of (12. 6p into (12. 7p results in the second order linear 
equation 

^m" - T (t) ij' + Aa (t) fi = (2.13) 

with 

, X ct' , / N , , ,0 d f a' d'\ , , 

rt = 2c + Ad, a(t)=ab-cd + d^ + -{ , 2.14 

a 2 \a d J 

which must be solved subject to the initial data 

/i(0) = 0, /i' (0) = 2a (0) 7^ (2.15) 

in order to satisfy the initial condition for the corresponding Green function; see the asymptotic 
formula (12.241) below. We shall refer to equation (12.131) as the characteristic equation and its 
solution fi{t) , subject to (12.151) . as the characteristic function. As the special case (12.131) contains 
the generalized equation of hypergeometric type, whose solutions are studied in detail in ^29j; see 
also [1], [28], [38], and [30]. 

Thus, the Green function (fundamental solution or propagator) is explicitly given in terms of the 
characteristic function 

^ ^ Q^^^ ^ ^i{a{t)x^+l3{t)^y+-i{t)y^+S{t)x+e{t)y+^,{t)) ^ (2.16) 

\/2T:%y. it) 



Here 



^^"4a(t)/x(t) 2a (t)' ^ ^ 

P it) = exp ^ (c (r) - 2d (r)) dr^ , (2.18) 

7W = u}^]^. f {c{r)-2d{r)) dr\ (2.19) 



/i(t)/i' (t) 



h (a'(t)) 

"t \ rt 

p-sm I — / 

ft it) 



6{t) = -l-^exp{- I {c{t) -2d{T)) drj J exp { j (c (A) - 2rf (A)) c/A ) (2.20) 



X ((/(-)- (r)) (r) + (r) ) dr, 



a (r) J 2a (r) 

e{t) = -^5(t) exp(^-£(c(r)-2rf(r)) rfr) (2.21) 
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' air) air) 



with 



HO) = {^y s{0) = -SiO), KiO) = 0. (2.23) 

We have used integration by parts to resolve the singularities of the initial data. Then the corre- 
sponding asymptotic formula is 

as t — s> 0+. Notice that the first term on the right hand side is a familiar free particle propagator 
(cf. ([31D below). 



i^ = Hit)^, ^ix,t)l^, = %ix) (2.25) 



By the superposition principle, we obtain an explicit solution of the Cauchy initial value problem 

.dip 
'dt 

on the infinite interval — oo < a; < oo with the general quadratic Hamiltonian as in (12.11) in the 
form 

POO 

iPix,t) = / Gix,y,t) ^/j,{y) dy. (2.26) 



This yields the time evolution operator (11.21) explicitly as an integral operator. Properties of similar 
oscillatory integrals are discussed in [36j . 

3. Some Special Cases 

Now let us consider several elementary solutions of the characteristic equation fl2.13p : more com- 
plicated cases may include special functions, like Bessel, hypergeometric or elliptic functions [T], 
[29j . [31J . and [ID]. Among important elementary cases of our general expressions for the Green 
function fl2:T6|) - fl2^ are the following: 

When a = fi/2m, b = c = d = f = g = 0, and fi" = 0, /i = tH/m, one gets 

/^-(^ I • 

which is the free particle propagator pT] . 
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For a particle in a constant external field, where a = H/2'm, b = c = d = g = and / = F/H 
constant, /i = th/m, the propagator is given by 



This case was studied in detail in [2], [0], [H], [U], and [53] ■ Here we have corrected a typo in 
[TT] ; see [37] for a complete list of known errata in the Feynman and Hibbs book. 



The simple harmonic oscillator with a = fl/2m, b = muj'^/2fi, c = d = f = g = and fi"+uj'^fi = 0, 
H = {h/muj) sin cut has the familiar propagator of the form 

G (x, y,t) = \ exp — (ix + v ) cosut — 2xy) , (3.3) 

^ ^ y 27rihsmujt ^\2hsmujt^^ ^' ^'j' ^ ' 



which is studied in detail at [4j, [13], [IS], [22], [21], [33- For an extension to the case of the forced 
harmonic oscillator including an extra velocity- dependent term and a time-dependent frequency, 
see [7], [8], [11] and [21]. 

Furthermore, an exact solution of the n-dimensional time-dependent Schrodinger equation for 
certain modified oscillator is found in [23j. In the one-dimensional case we get functions 

1 1 
a = -(l + cos2t), 6 = - (1 - cos2t) , c = 2rf = sin2t (3.4) 

and our characteristic equation (12.131) takes the form 

/i" + 2tant ^'-2^ = 0, (3.5) 

whose elementary solution is 

/i = cost sinht + sin t cosh t, (3.6) 
which satisfies the initial conditions (12.151) . Further, the corresponding propagator is given by 

G (x, y, t) = ^ = (3.7) 

^y 2-Ki (cos t sinh t + sin t cosh t) 

( [x^ — y"^) sin t sinh t -|- 2xy — (x^ + y^) cost cosh 

X exp — 

\ 2i (cos t sinh t + sin t cosh t) 

which was found in [23] as the special case n = 1 of a general n-dimensional expansion of the 
Green function in hyperspherical harmonics. We have showed that (13.71) is a generalization of the 
propagator for the simple harmonic oscillator; see Ref. [23] for more details. 



4. A Particular Solution of the Nonlinear Schrodinger equation 

The method of solving the equation (12.11) can be extended to the nonlinear Schrodinger equation 
with a general quadratic Hamiltonian of the form 

= W S + ^ - ' + ^ ^) (4.1) 

-/ it) X^ +ig{t)^ + h (t) ij, S>0, 
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where a{t) , b {t) , c{t) , d{t) , f (t) , g (t) , and h {t) are certain functions of time t only. A subs- 
titution of 

^ = ^ (x, t) = ! e^^(^'^'*\ = constant, (4.2) 

V/i (t) 

where S = S {x,y,t) = a (t) + jS (t) + 7 {t) y"^ + 6 (t) x + e (t) y + n (t) with the same relation 
(12. 6p . results in a modified system (I2.7I) - (I2.12I) below: 

— + b(t) + 2c (t) a + Aa (t) o? = 0, (4.3) 
at 

dd 

^ + (c(t)+4a(t)a(t))/3 = 0, (4.4) 

^ + a (t) {t) = 0, (4.5) 

^ + (c (t) + 4a (t) a (t)) 6 = fit) + 2a (t) g (t) , (4.6) 
dF 

- = {g{t)-2ait)6{t))f3{t), (4.7) 

^ = ^(t)5(t)-a(t)5^(t)-A^, (4.8) 

where only the last equation has an extra term which corresponds to the nonlinear term in the 
original Schrodinger equation (14. ip . Therefore equations ( 12.13p -( !2.16p solve once again the Riccati 
equation (14. 3p . however; in this case, we would like to use nonsingular initial conditions 

/i (0)7^0, /i'(0) = 2(2a(0)a(0) + rf(0))/i(0). (4.9) 

With these conditions, our system (I4.3p - (l4.8p can be solved again in terms of the characteristic 
function fj, (t) , thus giving us a particular solution of the nonlinear Schrodinger equation (14. ip in 
the form (14. 2p . corresponding to the initial data 

ijoix) = ^{x,t)\,^, = -^L= e<«(o)x^+/3(o).y+7(o)j/^+5(o).+e(o),+.(o))_ ^4_^Q) 

V/^(0) 

The details here are left to the reader. 
In the simplest case, 

dip T d^ijj 2 

i-— = —-;r + h\ih\^ih /i = constant, s > 0, (4.11) 

ot 2 ox^ 

the solution of the characteristic equation fi" = is fi (t) = fiQ + t/i^, /io > and the coefficients of 
the quadratic form are given by 

^(^) = vr^Tr-y Hi) = ^^^ 5{t) = ^^, (4.12) 

2 (/io + t/ii) /^o + tf^i /io + if^i 

l{i)=lo-7u — , ^ V £{t)=£o — — , 4.13 

K{t) = Ko- ^,^'fl ^ --L{t) (4.14) 
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with 



(l-s) 



{{fiQ + t^^f ' - /io , when s 7^ 1, 



ln(l + ^ 

^0 



when s = 1. 



Now, the hmiting case where /i^ — with /Iq > is given by 

Then \ip\ = is bounded at all times. Yet, when /i^ 7^ 0, one gets 

|^(x,t)| = ^^=, t>0, 

which is bounded if /Xq > and fi^ > 0, and blows up at a finite time to = ~/^o//^i A^i 
The same method shows that the Cauchy initial value problem 

.dtp l&^ip 



+ 



dt 2 9x2 

= Se{x-y) 

has the classical solution of the form 

i> = Ge {x,y,t) 

where 



hltpl ^ ijj, h = constant, s > 0, 

2' 



1 



Hie 



exp 



i {x - yY 



^Im (t + e) 



exp 



1 - s 



In ( 1 + - 



2e 

i{x — yY ih 
2{t + e) ~ 2^ 

-, when s 7^ 1, 
when s = 1 



£ > 



Xs (t) 



with (0) = 0. Here 

^L=o = '^e {x, y, 0) = 6s{x-y)^6{x- y) 
as £ — s> 0"^ in the distributional sense 



lim 

£^0+ 



Further details are left to the reader. 



Ge{x,y,^)^{y) dy = ^{x) 



In a similar fashion, the nonlinear Schrodinger equation 

M 2. , • 2. 2 , • • . +(0 

I—- = — cos t — — - + sm t X ip — t smr cost 2a;—— 

at ox^ \ ox 

(cf. |j23j) has a particular solution of the form 



ip = ip {x, t) 
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where n {t) = cos t cosh t + sin t sinh t and the coefficients are given by 

cos t sinh t — sin t cosh i z.^^x 

a (t) = — — , (4.25) 

2 (cos t cosh t + sm t smh t) 

P it) = r^—. r-^, (4.26) 

cos t cosh t + sm t smh t 

cos t sinh t + sin t cosh t (4 27) 

2 (cos t cosh t + sin t sinh t) ' 



(costcosht + sin t sinh t)"*^ — 1 



when s 7^ 1, 



(t) = <( 1 - s (4.28) 

— In (cos t cosh t + sin t sinh t) , when s = 1. 

Thus, initial function is the standing wave 

^|J,{x) = ^P{x,t)\,^, = e'^y. (4.29) 

These exact sohitions may be of interest in a general treatment of the nonlinear Schrodinger equation 
(see [17j|, yjj, ^26j, ^34j, [35], [H] and references therein). 

5. Motion in Perpendicular Magnetic and Electric Fields 

A particle with a spin s also has a magnetic momentum fj, with the operator 

p. = lj.s/s, (5.1) 

where s is the spin operator and /i is a constant characterizing the particle. For the motion of a 
charged particle in a uniform magnetic field H and an electric field E, both which are perpendicular 
to each other, the three-dimensional time-dependent Schrodinger equation 

ih— = H"^ (5.2) 

dt ^ ' 

has a Hamiltonian operator as in [20], namely. 



^ 1 f ^ eH \ ^ 1 1 

where p = —ifi^ is the linear momentum operator. The corresponding vector and scalar potentials 
are defined up to the gauge transformation. We use the original choice [20] for the vector potential 
A = —yH it) Bx and add a linear scalar potential ip = —{F{t)/e)y (see also [2T]). Then the 
uniform magnetic field H and the corresponding perpendicular electric field E are given by 

IdA _ydH (t) F (t) 

c dt c dt ^ e 

see llHl for more details. 



H = curl A = H{t) e„ E = - V<^ - -— = ^_lie^ + Sy] (5.4) 



Since (15.31) does not contain the other components of the spin, the operator commutes with 
the Hamiltonian H and the 2;-component of the spin is conserved. Thus the operator can be 
replaced by its eigenvalue = a in the Hamiltonian (15.31) . so we have 
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with a = —s, — ... , s—1, s. Then the spin dependence of the wave function becomes insignificant 
and the wave function in the Schrodinger equation fl5.2p can be taken as an ordinary coordinate 
function \1/ = \1/ fr, t, a) . 



It should be noted that, the Hamiltonian fl5.5l) does not contain the coordinates x and z exphcitly. 
Therefore the operators % and % also commute with the Hamiltonian and the x and z components 
of the linear momentum are conserved. The corresponding eigenvalues and Pz take on all real 
values from —oo to oo. In this Letter we consider the case when the magnetic field H and the 
electric force F are arbitrary functions of time t only. Using the substitution 

^ (r, t) = e*(^P-+^P--^^W)/« (y, t) (5.6) 

with 

dSp pI /id cp^ F{t) 

results in the one-dimensional Schrodinger equation of the harmonic oscillator driven by an external 
force in the ^/-direction 

with time-dependent quantities 



mc en [t) \lmujH\t) 

When the magnetic field is a constant H' = 0, the solution is well known (see [20], [21]). In the 
absence of external electric force F = 0, there exist discrete energy values corresponding to motion 
in a plane perpendicular to the magnetic field, namely, the Landau levels. See [20] , [21] and section 7 
for more details. 

In general, the following substitution 

i;iy,t)=xiv,t), V= ^~ ^° f ^ (5.10) 

an [t) 

gives the time-dependent Schrodinger equation of the form (12. ip as follows, 
where 

an {t)F{t) y'^{t) a'^ (t) 

which can be solved by the method discussed in the previous sections. In this case 

r (t) = + 2^ = (In {unal))' = 0, a (t) = \ul (t) 

and the corresponding characteristic equation (I2.13P coincides with the equation of motion for the 
classical oscillator with a time- dependent frequency 

/i^ + u;^(t)/XH = 0, /iH(0) = 0, /x^ (0) = (0) . (5.13) 
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The fundamental solution of (15.111) is given by 



where 



«H it) = (In;.^ it)) = ^al it) I (In/.^ (t)) 



^^^^ fiH{t)aH{0) finit)^ ujH{t) fi^it)^ H{t)' 

2 \fiHit)f^Hit) Jo V/^hM/ / 

(«) = f f H {r) F (t) dT, 
Jo 

|e| eh [t) mc Jq 
^' Jo \ f^H (t) /i^ (t) J {r) 



(t) = ^ (4°^ (t) + pJP (t) + pIkP (t)) , 

Jo 



2 ' ^' I^H{t)l^'Hit) Jo V/^^M 

F (r) (r) 



(ir. 
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.g.,)^ ^"'"i^-f'M'^^^^X.r (5.27, 



clS db special case of fl2.16l) - fl2.22l) . 



6. The Propagator in Three Dimensions 

By the superposition principle, the fundamental solution \E' = G (r, r', t) of the time-dependent 
Schrodinger equation of a particle with a spin in a uniform magnetic field 

ih^ - H] G{r,r',t) =0 (6.1) 



dt 

is given by the double Fourier integral of the form 



oo 



G(rr't) = ^ I I ^i{(x-x')p^+{z-z')p^-SF{t,p^,p^))/n 2) 

^ ' ' ^ (27r^)'aH(0) ' ' ^ ■ ^ 



a/f [t) an (0) 



Here we obtain from (15.71) that 



S,(t,,^,,^)^A.t,^[lS^^,r-^^[H(r) ,r. (6.3) 

By (15.141) the corresponding Green function is represented as 

Gh iv, V\ t) = , ^ e^^-^"'"' (6.4) 
^/2mfiH (t) 

with 

Sh iv, V', t) = OLE (t) rf + 13 H (t) W + 7^ (t) {rj'f + 5h (t) r] + Eh (t) v' + (t) (6.5) 

and 

^ y-yo jt) , ^ y' - yo (0) 

where ?/o is a linear function of [see (15. 9p for a definition of functions y^ (t) and a^j (t)]. Next, we 
are given 

hm G (r, r', t) = 5 (r - r') = 5 (x - a;') 5 (y - y') 5 {z - z') (6.6) 
as our initial data by the asymptotic relation (12.241) and the integral representation 

as the Dirac delta function. 

The quadratic in (16. 5p is also a quadratic polynomial in given as: 

Sh iv, V', t) = (y, y\ t) + p. S^^^ (y, y\ t) + (t) (6.8) 
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with coefficients 

^ ^ ' e^al (t) m (t) ^ e^an (t) an (0) H (t) H (0) e^aj, (0) (0) ^ ' 

, it) ^ ceg jt) ^ K^S it) 



nefiH{t)H{t) meajj [0) H {0) 2^m ' 

^H^iy^y'^t) = (^eal, it) Hit) + ean {tfal (0) H (0) + ^1^^) ^ ^^'^^^ 

/ (t) 2c7^ (^) gj;^ (t) ^ , 

\^eaj^ (t) Oh (0) H (t) eajj (0) (0) ma| (0) y ^ 

, c6P (t) ^ csf jt) ^ (t) 



and 



Hefiff (t) H (t) mea\ (0) H (0) 2^m 



(^H {t) clh [t) an (0) ajj (0) 

, SPit)^ , ^, , .Pit) 



(6.11) 



See equations fl5.9p and fl5.15l) - fl5.27l) for notation explanation. 

Hence, the double Fourier integral (16.21) can be evaluated in terms of elementary functions. So, 
integration over gives the free particle propagator of motion in the direction of magnetic field in 
the following fashion. 



oo 



m I im [z — z 

exp 



\2 



2Txint \ 2nt 

using the familiar elementary integral 

,<-^'+^'^)dz = \l-e-''''^- (6.13) 



oo 

e" 




see Refs. [3] and [30] . 

Similarly, integration over yields 



Gir,r',t) = ^^^^^^^^e.J'J^ f Hir) dr] (6.14) 

27THaHiO)^2fiHit)SPit) ^0 
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^ ^^U^ (y, y', t)) ' - 4sl?) (t) {y, y', t) ^ 



X exp 



V 



MSP it) 
2z/i4°) (t) V'' ~ " e io H (r) 



with 

H I he 



an (0) 



muHiO) Y |e|//(0) 
Moreover, the quadratic form in this expression for the 3D-propagator 

Q (y, y\ t) = (Sg) {y, y\ t))' - ASP (t) SP (y, y', t) 
(or the discriminant of (16.81) ) can be rewritten as 

Q (y, y', t)=A (t) y^ + B {t) yy' + C {t) {y'f + D {t)y + E {t)y' + L (t) , 



where 



2 dy"^ ' dydy' ' 2 ' 

dQ 

r, = 

dy' 



D 



y=y'=0 



L=Q\ 



y=y'=o 



y=y'=0 ' 



Our final result is 

A{t) 



dy 

c2 {(3j, (t) - Aan (t) 7^ (t)) 2c/5^ (t) (t) 



+ 



e'^ajj {t) al {0) H^O) hefiH (t) an {t) an (0) H {0) 



B{t) 



1t?^i]j (t) meajj (t) aj^ (0) H (0) /ima|^ (t) 
2c' {Pj, (t) - Aau (t) iH it)) 2cPh (t) 4^ (t) 



e^al^ (t) (0) if^ (q) ^e/i^ (t) (t) a^^ (0) H (0) 



meal (t) (^h (0) H (t) hcfiH (t) ajj (0) H (0) ^m/x^ (t) (0) 
2c/5^ (t) eg) (t) 2/5^ (t) (t) 



mean (t) af^ (0) H (0) HmaH (t) an (0) ' 
_ c^(/3^(t)-4a^(t)7^(t)) , 2cPHit)e^^{t) 



e'^ajjit) al {0) H^t) mean (t) a^iO) H {t) 
4^ (^)) 4c7h (i) it) 47h (t) (t) 



m^af^ (0) ^e/i^ (t) ajj (0) if (t) 2/ima|^ (0) 
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^ ' me^ajj (t) ajj (0) H (t) H (0) me^an (t) a% (0) (0) ^ ^ 



HmefiH (t) af^ (0) H (0) /imea|^ (t) H (t) 



hmeaH (t) au (0) H (0) h^mnu (t) 

2c' Ph (t) (t) Ac'^H (t) 5P (t) 



he^^H (t) (t) an (0) H (t) H (0) He^HH (t) < (0) (0) 
2c6P {t)6P (t) 



n^ei^l (t) H (t) ' 



p.. ^ 2c^/3^ (t) SP it) 4c^^ (t) 4°) it) 

He^f^H (t) an (t) an (0) ^ (t) He^f^H (t) a% (0) H (t) H (0) ^ ^ ^ 

2c [SP it) eP it) - 26P it) eP {t)) 



+ 



Hmefifj (t) (0) if (t) Hmean {t) an (0) if (t) 

2c7h (t) 4'^ (t) ^ 4^ (^) 4^ (^) - {t) 4^ (t) 



Hmeajj (0) H (0) ^m2a|^ (0) 

Ac^aH(t)ePit) 2c'(3H{t)eP{t) 



' me^af^ (t) a]j (0) ^ (t) me^aH (t) a% (0) H (t) H (0) 

2c6P (t) 6P (t) 



ea% (0) H (0) ' 



^ {SP it)) c'{ePit)) {4^ (t)^ 

, 2c'6Pit)ePit) ^ c6P(t)4'(t) , c4°^(t)K«(t) 



" hme^fiH (t) ajj (0) H (t) H (0) /i^^e/x^ (t) H (t) nw?ea\j (0) H (0) 

2c2aff (t) /cg^ it) 2c2/j^ (t) /eg) it) 2c'^H it) 4^ it) 

Hme^ajj (t) ii2 (t) Hme^an it) an (0) H (t) H (0) te^a^^ (0) H'^ (0) 

2c4) (t) Kg) (t) _ 2c5i^) (t) it) _ 4^ it) 4^ jt) 
n^mei^H it) H (t) hmWHiO)H{0) f?w? 

The details are left to the reader. 

Finally, by the superposition principle, a general solution of the Cauchy initial value problem in 
i?^ subject to the initial data 

* (r-, t) |,=o = * {r, 0) = (x, y, z) (6.24) 
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has the form 

^(r,t)= / G{r,r',t) ^(r',0) dx'dy'dz. (6.25) 

This gives us the time evolution operator (11.21) exphcitly for a motion of a charged particle in a 
uniform magnetic field and also a perpendicular electric field with a given projection of the spin 
Sz = (J m the direction of magnetic field. 



7. Two Examples 



7.1. Motion in a Constant Magnetic Field. The simplest case occurs as if' = F = 0, which 
implies //^ + w^fifj = with fi^j = sin {unt) . Thus, 

(^nit) =-fH (t) = 1 cot (unt), (3H{t) = — ■ / (7.1) 
2 sm {uJut) 



m _ _ .(0) _ .(1) _ ^(0) _ ^(1) _ ^(2) _ r. 

Up — U — tp — tjj- — h,p — h,p — h,p[ — U, 



and 



52>fe..',*)^-i^ta„(^) (7.2) 
S2' (!/,!/'.«) = ^ (cot {^„t)y'- + cot (u^Ht) {yf 



where the discriminant is 



Q {y, y') = (5« {y, y', t)) ' - 4^°) (t) {y, y', t) = ^^{y-y'f. (7.3) 
Hence, the Green function is given by 



X exp 



Ah 

See pi| for more details. 



{{x - x'f + (y - y'f) cot {uHt/2) - 2^ - x') {y + y') 



7.2. A Linear Magnetic Field. Now consider the case H {t) = Hq + tHi, where Hq and Hi are 
constants. The characteristic equation (I5.13P becomes a special case of the Lommel equation [1], 
[28], [38], which can be solved in terms of Bessel functions [40] of orders u = ±1/4. It then follows 
that 



(^-1/4 (■ 



\e\H'o\j f \e\H^t) \ f \e\ \ ^ f \e\ H^t) 



2mcHiJ ' \ 2mcHi J ' \2mcHiJ ' \ 2mcHi 
with 

.2 /u^u \ 3/2 

(7.6) 



dt ~ w?c^Hi V 2 j 
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X 



J (\e\Hl\ ( \e\H^t) \ (\e\Hl\ ( \e\H^t) 



A general expression of the propagator is given above by (16.141) . In the case when F = 0, one can 
simplify = e 
cases elsewhere. 



simplify 5^"* = e:^^ = np = = 0. We shall elaborate on these and some other interesting special 
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